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A calculation of the leading isospin-breaking contributions to the B — > K*'y decay 



■ amplitudes based on the QCD factorization approach is presented. They arise at 



order A/mb in the heavy-quark expansion and are due to annihilation contribu- 
tions from 4-quark operators, the chromo-magnetic dipole operator, and charm 
penguins. In the Standard Model the decay rate for — > K*^^ is predicted to 
be about 10-20% larger than that for B~ — > K*~'y. Isospin-breaking effects are 
a sensitive probe of the penguin sector of the effective weak Hamiltonian. New 
Physics models in which the hierarchy of B ^ K*j decay rates is either flipped or 
greatly enhanced could be ruled out with more precise data. 



1 Introduction 



The study of radiative decays based on the flavor-changing neutral current transition 
6 — s> 57 is of crucial importance for testing the flavor sector of the Standard Model 
and probing for New Physics. Whereas the inclusive mode B Xg'y can be analyzed 
using the operator product expansion, it is usually argued that exclusive decays such 
as B —>■ K*'y and B ^ p'j do not admit a clean theoretical analysis because of their 
sensitivity to hadronic physics. However, it has recently been shown that in the heavy- 
quark limit the decay amplitudes for these processes can be calculated in a model- 
independent way using a QCD factorization approach [|l], 0, ^, which is similar to the 
scheme developed earlier for non-leptonic two-body decays of B mesons 

To leading order in A/rrih (and neglecting isospin violation in the B —>■ K* form 
factors) one finds that the amplitudes for the decays B^ K*^'y and B^ — > K*^'y 
coincide. Spectator-dependent effects enter at subleading order in the heavy-quark ex- 
pansion. In this Letter the QCD factorization approach is used to estimate the leading 
isospin-breaking effects for the B K*'~f decay amplitudes, the most important of which 
can be calculated in a model-independent way. 

Experimental measurements of exclusive B — > K*'~f branching ratios have been re- 
ported by the CLEO, Belle and BaBar Collaborations, with the results (averaged over 
CP-conjugate modes): 



4.55l°:6l±0.34 [5] 

10^ Br(B" ^ K*°7) = { 4.96 ± 0.67 ± 0.45 || 

4.23 ±0.40 ±0.22 0] 

3.76i°:i±0.28 i 

10' Br(5- ^ K*-7) = <{ 3.89 ±0.93 ±0.41 § 

3.83 ±0.62 ±0.22 M 



(1) 



The average branching ratios for the two modes are Br(i?° — > K*'^'y) = (4.44±0.35) ■ 10~' 
and Br(_B~ K*~'-)) = (3.82 ± 0.47) ■ 10~^. When corrected for the difference in the 
5-meson lifetimes, tb- /t^o = 1.068 ± 0.016 these results imply 

Although there is no significant deviation of this quantity from zero, the fact that all 
three experiments see a tendency for a larger neutral decay rate raises the question 
whether the Standard Model could account for isospin-breaking effects of order 10% in 
the decay amplitudes. 
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2 Isospin-Breaking Contributions 



In the Standard Model the effective weak Hamiltonian for 6 — 57 transitions is 

1=3,. ..,8 



V ^ p=u,c ^ i=3 8 ^ 



where A^**^ = ^^H^fc products of CKM matrix elements, Q12 are the current-current 
operators arising from W exchange, Qs^...^ are local 4-quark penguin operators, and 
Qj and Qs are the electro-magnetic and chromo-magnetic dipole operators. (We adopt 
the conventions of [Q; in particular, Ci ~ 1 is the largest coefficient.) The Wilson 
coefficients Ci and the matrix elements of the renormalized operators Qi depend on the 
renormalization scale fi. 

At leading power in A/m^, and neglecting the tiny contribution proportional to A^'*^ 
the B K*'y decay amplitude is given by 

^^lead = ^ Xi'^a', {K*{k,r^h{q,e)\Qr\B) , (4) 

where the next-to-leading order (NLO) result for the coefficient = C7 + . . . is given in 
Eq. (42) of 0]. At this order contributions to the quantity Aq- arise only from isospin 
violation in the B K* form factors. (The phase-space difference between the two 
decays is an effect of order (A/m^)^ and contributes the negligible amount —4 x 10~^ 
to Aq-.) To first order we find 

A soft _ . _ Tf^^ _ mg - rriu / T\ ^\ _ ^ , . 

where Tf^^* is a form factor in the decomposition of the B K* matrix element of 
the tensor current S(T^j,(1 + 75)6 evaluated at zero momentum transfer. We have used the 
form-factor predictions of and assumed that the ratio of isospin to ?7-spin violation 
for the transition form factors scales approximately as the corresponding ratio of light 
quark masses. Although this estimate is rather uncertain, we believe it indicates that 
the "soft" isospin-breaking effect is negligible and could not account for a value of Aq- 
as large as 10%. 

At sub leading order in A/rrif, "hard" isospin- violating effects appear in the form of 
spectator-dependent interactions. The leading contributions arise from the diagrams 
shown in Figure For a first estimate of these effects we adopt a simplification of 
the usual NLO counting scheme, in which we neglect terms of order as Cs^,,,^ while 
retaining terms of order a^Ci^g- This is justified, because the penguin coefficients C^^,,,^ 
are numerically very small. Also, it is a safe approximation to neglect terms of order 
as /^i^-' given that lAi'V^c'^l ~ 0.01-0.02 is very small. It then suffices to evaluate 
the contributions of the 4-quark operators shown in the first diagram at tree level. The 
terms neglected in this simplified NLO approximation will be estimated later. 
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Figure 1: Spectator-dependent contributions from local 4-quark operators (left), the 
chromo- magnetic dipole operator (center), and the charm penguin (right). Crosses 
denote alternative photon attachments. 



The isospin-breaking contributions to the decay amplitudes can be parameterized as 
Aq = hq A\ea.Ai wherc Q is the flavor of the spectator antiquark in the B meson. To leading 
order in small quantities Aq- is then given by 



Ao- = Reipd - hu) . 



(6) 



The QCD factorization approach gives an expression for the coefficients hq in terms of 
convolutions of hard-scattering kernels with light-cone distribution amplitudes for the 
K* and B mesons. When light quark masses are neglected, the leading and subleading 
projections for a transversely polarized vector meson with momentum k and polarization 
r/ are H O 



{K\k,r,)\-s{-z)^...q{z)p\Q) 
fK*mK 



)l3a 



+ 



Jo 

T]*- Z 



x] 



2 ^/lUpCT 



r/*^F^'^(7^5)/3a / dxe'^''-^g'f\ 
Jo 



x] 



+ 



k ■ z 



(v) j ^ 

91 



(7) 



where = 0, and the ellipses on the left-hand side indicate a string operator required 
to make the non-local matrix element gauge invariant. The variable x is the longitudinal 
momentum fraction of the strange quark, and C = (1 — x) — x. All four distribution 
functions are normalized to 1. The asymptotic form of the leading- twist amplitudes 
0_Ljj(x) is 'oxx with a; = (1 — x). In the approximation where 3-particle distribution 
amplitudes of the kaon are neglected, the functions 0|| and 5'^''*'' are related to each 
other by equations of motion [|T^, |rT|. Using these relations we find that (the prime 
denotes a derivative with respect to x) 



.r(.)-^ 



(8) 
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which can be used to ehminate 0|| from our results. This means that we neglect con- 
tributions to the quantity K2 arising from {gqs) Fock states of the kaon, which are 
typically found to be suppressed with respect to two-particle contributions of the same 
twist 0. This approximation is justified, because numerically the effect of K2 is about 
four times smaller than that of Ki. Including the 3-particle Fock states would, however, 
not invalidate factorization. 

The leading-twist projection onto the B meson involves two distribution amplitudes 
^Bi{0 ^iid $B2(0 0) 0! where ^ = 0{K/mij) is the light-cone momentum fraction of 
the spectator quark projected onto the direction of the kaon. The first inverse moment 
of the function defines a hadronic parameter Xb via 



i 



TUB 



(9) 



The function ^b2{0 does not enter our results. Note that 3-particle Fock states of the 
B meson only contribute to the quantity Ki at order ^^(mf,) C^fii^rrih) |jl3[ and thus can 
be neglected in our approximation scheme. 

The leading spectator-dependent contributions summarized by the coefficients bq can 
be written as 



mbT{ 



rrib 



GXBmB 



where ruf, denotes the running 6-quark mass. The product rUbTf^^ is scale inde- 

imply 



pendent at NLO. The heavy-quark scaling laws Jb ~ Trijj and 
that bq scales like A/rUh, and hence the spectator-dependent corrections contribute at 
subleading order in the heavy-quark expansion. However, because of the large numerical 
factor 127r^ in the numerator the values of bq will turn out to be larger than anticipated 
inii- 

The dimensionless coefficients Ki are given by (A^ = 3 and Cp = 4/3 are color factors) 



■^B^K* 



-3/2 



Ca 



9l 

N 



+ ri 



C2 
N 



Ch\ Cf as 



(11) 




+ r2 



where ri and r2 are the residual NLO corrections neglected in our approximation scheme. 
The quantities 



G , (Sr) 







1 0±(x) 
dx——, 

ax — G(Sc, X) 



3x 
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Xi. = J^'dxM^)^^ (12) 

with Sc = {nic/mjjY are convolution integrals of hard-scattering kernels with the meson 
distribution amplitudes, and 

G'(s,x) = — 4 / du uu \\i{s — uux — ie) (13) 
Jo 

is the penguin function. The terms in K2 arising from the current-current operators Q12 
only contribute for g = m, as indicated by the symbol 5qu- Since no distinction between 
CP-conjugate modes is made in the experimental determination of Aq- we only need 
the real part of the ratio of CKM parameters, given in terms of Wolfenstein parameters 
as Re(A(fV^i'^) = = 0.011 ± 0.005. 

The first three convolution integrals above exist for any reasonable choice of the 
distribution amplitudes. This shows that, to the order we are working, the QCD fac- 
torization approach holds at subleading power for the matrix elements of the 4-quark 
and current-current operators (including penguin contractions). Factorization of these 
matrix elements is indeed expected to hold to all orders in perturbation theory as a 
consequence of color transparency |l[]. This implies, in particular, that long-distance 
contributions to the first diagram in Figure |1|, which have been analysed using QCD 



sum rules in [|T^, |T5|], must be suppressed by at least two powers of A/mt and so do 
not contribute to the order we are working. On the other hand, if the function <i)±{x) 
vanishes only linearly at the endpoints (as indicated by its asymptotic behavior), then 
the convolution integral X±_ suffers from a logarithmic endpoint singularity as x — > 1, 
corresponding to the region where the light spectator in the K* meson is a soft quark. 
This indicates that at subleading power factorization breaks down for the matrix element 
of the chromo-magnetic dipole operator Qg (second diagram in Figure |I|). In the phe- 
nomenological analysis below we regulate this singularity by introducing a cutoff such 
that X < 1 — Kh/niB, where Kh = 0.5 GeV is a typical hadronic scale. Since X± is 
dominated by soft physics, we assign a large uncertainty to this estimate. 

The Wilson coefficients C3 ...^g entering the expressions for Ki and K2 must be evalu- 
ated at NLO, while the remaining coefficients can be taken at leading order. Throughout, 
the two- loop expression for the running coupling ^^(/i) is used. The explicitly scale- 
dependent terms in the expressions for the coefficients Ki arise from the charm-penguin 
diagrams. They are necessary to cancel the scale and scheme dependence of the Wilson 
coefficients C3,...,6- In a complete NLO calculation there would be additional logarithmic 
terms proportional to C*3,...,6; which can be deduced using the renormalization group 
(RG). This gives for the remainders 



ri 



r2 



^ C3 + ^ ( C4 + Ce ) - 8 ( iVCe + C5 



1^0 



44 4 

-^C,--nf{C, + C, 



In^-F..., (14) 

1^0 
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where nf = 5, and fiQ = 0{mh) is an arbitrary normalization point. The sensitivity to /iq 
provides an estimate of the residual non-logarithmic NLO terms denoted by the ellipses, 
whose calculation is left for future work. After the addition of the pieces the quantity 
K2 is RG invariant at NLO. In the case of Ki a scale dependence remains, which cancels 
against the scale dependence of the tensor decay constant, f^'if^) ~ [c(s{fJ^)]^^^^° , and 
of the running 6-quark mass. 

An important element of the phenomenological analysis are the convolution inte- 
grals. We adopt the shapes of the light-cone distribution amplitudes obtained (at sec- 
ond order in the Gegenbauer expansion) using QCD sum rules and vary the ampli- 



tude parameters within their respective error ranges. This leads to F± = 1.21 ± 0.06, 
Gx(sc) = (2.82 ± 0.20) + (0.81 ± 0.23)i and H^{sc) = (2.32 ± 0.16) + (0.50 ± 0.18)i, 
where mc/mh = 0.26 ± 0.03 has been used for the ratio of quark masses. Next, we find 
X± = (3.44 ± 0.47) X - (3.91 ± 1.08), where X = \n{mB/Ah) (1 + ge''^) parameterizes 
the logarithmically divergent integral Jq dx/x. Following ^ we allow g < 1 and an ar- 
bitrary strong-interaction phase ip to account for the theoretical uncertainty due to soft 
rescattering in higher orders. The above results for the convolution integrals refer to a 
renormalization point of v^GeV. (The expressions for given earlier are valid if the 
convolution integrals are normalized at a fixed scale.) 

Further input parameters are mf,(mb) = 4.2 GeV, the decay constants /b = (200 ± 
20) MeV 0, ifg, = (226 ±28) MeV and = (175 ±9) MeV (at fx = y/bGeV) 



and the parameter = (350 ± 150) MeV defined in terms of the first inverse moment of 
the 5-meson distribution amplitude $^1 HI- A dominant uncertainty in the prediction 
for Aq- comes from the tensor form factor Tf^^* , recent estimates of which range from 
0.32 [|1|] to 0.38 ± 0.06 0. On the other hand, a fit to the B K*-f branching 
fractions yields the lower value 0.27 ±0.04 To good approximation the result for Aq- 
is inversely proportional to the value of the form factor. Below we take Tf~*^* = 0.3 
(at /i = mfe) as a reference value. 

From the diagrams in Figure |I] it is seen that in all cases the operators are probed 
at momentum scales of order fi ~ rrih. Hence, following common practice we vary the 
renormalization scale between mh/2 and 2mf,. The result for Aq- is shown in Figure ^ 
The width of the band reflects the sensitivity to input parameter variations. The three 
curves correspond to different choices of the scale /xq in the expressions ( [T^ ) for the 
remainders rj. The excellent stability under variation of both /i and /io shows that our 
approximation scheme captures the dominant terms at NLO. 

Combining all sources of uncertainty we obtain 

Ao- = (8.0j:^i)%x^^. (15) 

The three largest contributions to the error from input parameter variations are due to 
the divergent integral X± (±1.2%), and the decay constant fs (±0.8%). 
The perturbative uncertainty is about ±1%. Our result is in good agreement with the 
current central experimental value of Aq- including its sign, which is predicted unam- 
biguously. By far the most important source of isospin breaking is due to the 4-quark 
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Figure 2: Prediction for the quantity Aq- as a function of the renormahzation 
scale, assuming T^^^* = 0.3. The dark hues refer to /io = mi, (sohd), mb/2 (upper 
dashed) and 2mh (lower dashed). The band shows the theoretical uncertainty. 

penguin operator Qq, whose contribution to Aq- is about 9% (at /i = /iq = mb). The 
other terms are much smaller. In particular, the contribution of the chromo-magnetic 
dipole operator, for which factorization does not hold, is less than 1% in magnitude and 
therefore numerically insignificant. Hence, the most important isospin-breaking contri- 
butions can be calculated using QCD factorization. It follows from our result that these 
effects mainly test the magnitude and sign of the ratio Cq/Ci of penguin coefficients. 

3 New Physics 

Because of their relation to matrix elements of penguin operators, isospin-breaking effects 
in B ^ K*'y decays are sensitive probes of physics beyond the Standard Model. For 
example, scenarios in which the sign of Aq- is flipped could be ruled out with more 
precise data. In addition, in certain extensions of the Standard Model there exist local 
4-quark operators yielding an isospin-breaking contribution to the decay amplitudes at 
leading power in the heavy-quark expansion. Precise measurements of radiative decay 
rates would tightly constrain the corresponding Wilson coefficients. Here we confine 
ourselves to a brief illustration of the most interesting potential New Physics effects in 
B 7^*7 decays. A more detailed study will be presented elsewhere. 

New Physics effects arising at some high energy scale manifest themselves at low 
energy through new contributions to the effective weak Hamiltonian. As a first, popular 
example in which the operator basis is not enlarged, consider the minimal supersymmet- 
ric Standard Model (MSSM) with minimal flavor violation, and with contributions to 
the B Xs'j decay rates that are enhanced in the large-tan j3 limit taken into account 
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beyond leading order ^ |^. In this scenario new contributions to Q^^,,,^ and Qg 
are too small to have a significant effect. For low tan/3, 'Re{Ci) is negative as in the 
Standard Model. However, for large tan (3 the coefficient Re(C7) can take on both posi- 
tive or negative values. Positive values, which would ffip the sign of Aq-, become more 
probable as tan/3 increases (see, e.g., p2|). With more precise data, isospin breaking in 
B K*'~f decays could rule out significant regions of the MSSM parameter space at large 
tan/3. We note that the sign of Re(C7) can also be fiipped in supersymmetric models 
with non-minimal fiavor violation, independently of tan (3, via gluino-down-squark loop 
graphs. 

In a more general scenario. New Physics effects can be parameterized in terms of 
additional contributions to the coefficients of the dipole operators and Qs, contribu- 
tions involving dipole operators and Qs with opposite chirality than in the Standard 
Model, and a plethora of local 4-quark operators, the most general set of which can be 
taken as (a factor of ^ A^*-* is included here for convenience only) 

K'ff7p' = %AWE E + (16) 

q=u,d,... r=ri(g.r2 

where Op = sTibqT2q and Op = SiTibj qjT2qi. The Dirac structure F = Fi r2 can be 
any of the ten combinations (F ip A) O (l^ =f A), {VtA)0{V±A), (5 ± P) (g) (S* =f P), 
{S ±P) (g) {S ± P), and Tl,r (S) Tl,r. Here, as usual, S" = 1, P = 75, 1/ = 7^, A = 7^75, 
and Tl^r = (t^i/(1 =F 75)- As an example, note that such an extensive list of 4-quark 



operators can arise in general supersymmetry models [^, |J . 

In a generic New Physics model the Wilson coefficients Cp and Cp need not be fiavor 
independent [|^. Therefore, isospin-breaking effects can arise even if the photon in the 
first diagram in Figure |I| is emitted from the b or s-quark lines. Taking this possibility 
into account, we find that at leading order the contributions to the coefficients Ki are 

.NP _ (n I q \ f Qd ^ \n 



-^2 — ^(V-A)®(V~A) r)^{S-P)®(S+P)-> \^') 



where Cp = -|- Cp/A^, and we have introduced the new convolution integral 

rfx = 0.84 ± 0.06 . (18) 

3a; ^ ' 

The New Physics contribution from the chromo-magnetic dipole operator is included in 
(p!7|) despite its 0{as) suppression, because it could potentially be large in models with 
strongly enhanced Og. In analogy with (|TUp, we define a quantity b^^ with corresponding 
coefficients Kf^ given by equivalent expressions with all Wilson coefficients replaced by 
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their opposite-chirality counterparts. Its contribution to the decay amphtude is Ag = 
b^^ Aicad, where ^lead is defined in analogy with @ in terms of the matrix element of 
the opposite-chirality operator Qy. Because this amplitude does not interfere with the 
leading Standard Model amplitude for B K*'y its effect is likely to be suppressed. 

A remarkable fact is that there can exist local 4-quark operators yielding a leading 
(i.e., not power-suppressed) contribution to the decay amplitudes, as indicated by the 
factor ~ itlb/Xb in the first term in Kf^ and K^^ , which compensates the power 
suppression from the prefactor f^,/mb in (|T0|). The origin of these terms can readily be 
understood by noting that the tensor operators scr^,^(l±75)gg(T^'^(l±75)6 have a leading- 
twist projection onto the K* meson. When Fierz-transformed into the basis used above, 
such operators turn into tensor operators and operators with structure {S±P)®{S±P). 
The fact that such operators enter the amplitudes for exclusive radiative B decays with 
large coefficients means that future, precise measurements of radiative decay rates will 
provide tight constraints on the corresponding Wilson coefficients. 

4 Summary 

We have presented a model-independent analysis of the leading isospin-breaking contri- 
butions to the B K*'y decay amphtudes. In the Standard Model these effects appear 
first at order A/m^ in the heavy-quark expansion. They can be expressed in terms of 
convolutions of hard-scattering kernels and meson light-cone distribution amplitudes of 
leading and subleading twist. We have evaluated these contributions at NLO in pertur- 
bation theory, neglecting however some numerically suppressed 0{as) terms. With the 
exception of the matrix element of the chromo- magnetic dipole operator, whose contri- 
bution is numerically small, factorization of the leading isospin-breaking contributions 
is expected to hold to all orders in perturbation theory. 

To our knowledge, the analysis presented here provides the first example of a quan- 
titative test of QCD factorization at the level of power corrections. As such, it lends 
credibility to the idea of factorization as a leading term in a well-behaved expansion 
in inverse powers of the 6-quark mass. Our prediction for the magnitude and sign of 
isospin-breaking is in good agreement with the present central experimental value of this 
effect. If this agreement persists as the data become more precise, it will be possible to 
place novel constraints on flavor physics beyond the Standard Model. 
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